Abstract. Results of a systematic study of substructure 
in X-ray surface brightness distributions of a combined 
sample of 470 REFLEX+BCS clusters of galaxies are pre- 
sented. The fully automized morphology analysis is based 
on data of the 3rd processing of the ROSAT All-Sky sur- 
vey (RASS-3). After correction for several systematic ef- 
fects, 52 ± 7 percent of the REFLEX+BCS clusters are 
found to be substructured in metric apertures of 1 Mpc 
radius (Hq — 50 kins" 1 Mpc -1 ). Future simulations will 
show statistically which mass spectrum of major and mi- 
t-H j nor mergers contributes to this number. Another impor- 

■ tant result is the discovery of a substructure-density re- 
lation, analogous to the morphology-density relation for 

■ galaxies. Here, clusters with asymmetric or multi-modal 
X-ray surface brightness distributions are located pref- 

' erentially in regions with higher cluster number densi- 
ties. The substructure analyses techniques are used to 
compare the X-ray morphology of 53 clusters with ra- 
! dio halos and relics, and 22 cooling flow clusters with the 
• REFLEX+BCS reference sample. After careful equaliza- 
' tion of the different 'sensitivities' of the subsamples to 
C*~) . substructure detection it is found that the halo and relic 
| sample tends to show more often multi-modal and elon- 
£^ . gated X-ray surface brightness distributions compared to 
' the REFLEX+BCS reference sample. The cooling flow 
. clusters show more often circular symmetric and unimodal 
| distributions compared to the REFLEX+BCS and the 
O ■ halo/relic reference samples. Both findings further sup- 
port the idea that radio halos and relics are triggered by 
merger events, and that pre-existing cooling flows might 
be disrupted by recent major mergers. 

^ ' Key words: clusters: general - clusters: substructure 
• i— i . 
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1. Introduction 

Current structure formation scenarios suggest an hierar- 
chical growth of cosmic objects (e.g., Peebles 1980, White 
1996) where the merging of subclumps turns out to be 
a fundamental process. Linear theory predicts that clus- 
tering in critical density universes continues to grow up 
to present-day redshifts, and structure formation begins 
to decline already at z « il^ 1 — 1 in low density uni- 
verses, where f2 is the present value of the normalized 
cosmic matter density. Compared to the Einstein-de Sit- 
ter case, clusters in low-f2 universes are thus expected to 
be more relaxed, less substructured, and less elongated as 
shown by the simulations of the Virgo Consortium (see 
Thomas et al. 1998) and the simulations of Evrard et al. 
(1993), Crone, Evrard & Richstone (1996), Mohr et al. 
(1995), etc. The frequencies of subclumps (substructures) 
and elongations are thus useful statistical quantities with 
a direct relation to cosmology. See also the analytic work 
of, e.g., Richstone, Loeb & Turner (1992) and Lacey & 
Cole (1993). Note, however, that the effects are difficult 
to measure so that the resulting constraints on structure 
formation models are presently less stringent compared 
to, e.g., a direct measurement of the power spectrum of 
cluster number density fluctuations as presented in, e.g., 
Schuecker et al. (2001). 

Cluster mergers produce moderately supersonic 
shocks, compressing and heating the intracluster gas, and 
increasing pressure and entropy. This can be measured as 
local distortions of the spatial distribution of X-ray tem- 
perature and surface brightness (e.g., Schindler & Miillcr 
1993, Rocttigcr et al. 1997). Moreover, mergers affect clus- 
ter X-ray luminosity, magnetic field, and electron/ion non- 
equipartition (see, e.g., Schindler & Miillcr 1993, Roet- 
tiger, Stone & Burns 1999, Takizawa 2000, Ricker & 
Sarazin 2001). Therefore, studies of individual clusters 
provide a wealth of new and useful information on the 
physics of the merger process, e.g., from spatially resolved 
X-ray spectroscopy. 
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A more detailed study would thus ideally use deep 
pointed X-ray observations to get high signal-to-noise tem- 
perature and surface brightness maps. Unfortunately, X- 
ray data from public archives generally allow precise mea- 
surements on interesting but not necessarily representa- 
tive clusters of galaxies. For the majority of the clusters 
needed for a statistically representative sample, detailed 
studies are at present and in the near future not achiev- 
able. What is needed for cosmological investigations is a 
systematic and comprehensive study of substructure char- 
acteristics for a large set of clusters compiled in an homo- 
geneous way. 

Concerning the projected X-ray surface brightness dis- 
tribution, merger events cause multiple X-ray luminos- 
ity peaks, isophote twisting with centroid shifts, elonga- 
tions, and other irregularities. One of the first system- 
atic studies of such cluster X-ray morphologies was un- 
dertaken by Jones & Forman (1999) using the spatial sur- 
face brightness distributions of targeted and serendipitous 
clusters obtained with the Einstein imaging proportional 
counter (IPC). Their visual classifications of the X-ray iso- 
intensity contours of 208 clusters with z < 0.15, are useful, 
but subjective. For refined statistical analyses it should 
be supplemented by a more homogeneous sample selec- 
tion and by a more objective, i.e., quantitative method 
to analyse the cluster morphology. Along this line, Mohr 
et al. (1995) measured emission- weighted variations of the 
cluster centroid for substructure detection, but only for a 
comparatively small sample of 65 bright Einstein clusters. 
The remaining studies not mentioned here are restricted 
to even smaller cluster samples. 

One possibility towards a systematic study of substruc- 
ture is offered by the large area ROSAT All-Sky Survey 
(RASS, Triimper 1993, Voges et al. 1999). However, both 
the energy resolution of 40 percent FWHM at 1 keV and 
the limited energy range 0.1 < E < 2.4 keV of the RO- 
SAT position-sensitive proportional counter (PSPC) limit 
temperature measurements of the intracluster gas, espe- 
cially for kT > 4keV. This and the fact that the small 
number of RASS-3 photons counted in the direction of a 
typical galaxy cluster in our samples do neither allow the 
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determination of useful temperature maps nor the proper 
quantification of the amount of substructure, e.g., in the 
form of mass estimates of the individual subclumps for 
the majority of a fair cluster sample (but see Sect.||). The 
present analysis thus gives more weight to a representa- 
tive study of a large sample of clusters, mainly restricted 
to substructure detection and significance determination 
in terms of significance for the degree of substructure for 
a whole sample or for specific subsamples. 

In a series of papers we want to study the morphol- 
ogy of a large set of X-ray clusters in a systematic way. 
The present paper investigates the morphology of galaxy 
clusters using the X-ray surface brightness distributions 
extracted from data obtained in the course of the 3rd pro- 
cessing of the ROSAT All-Sky Survey (RASS-3). Further 
papers will concentrate on the study of alignment effects of 
the major cluster axes, and on detailed comparisons with 
large scale, high resolution numerical simulations. The 
morphological analyses are mainly restricted to galaxy 
clusters of the two largest and almost complete X-ray 
cluster surveys finished to date: the ROSAT-ESO Flux- 
Limited X-ray (REFLEX) cluster survey (Bohringer et al. 
2001) in the South, and a survey yielding the Brightest 
Cluster Sample (BCS, Ebeling et al. 1998) in the North. 

In the present investigation we also describe more de- 
tailed studies of specific cluster types with so-called ra- 
dio halos and radio relics, and signatures of cooling flows. 
Reviews concerning the latter cluster type can be found 
in, e.g., Fabian, Nulsen & Canizares (1984), and Fabian 
(1994), but see also Makishima et al. (2001). A few com- 
ments concerning the former two cluster types will be 
given in the following (see also the conference proceedings 
edited by Bohringer, Feretti & Schuecker 1999). 

Willson (1970) discovered diffuse, cluster-wide, steep 
spectrum synchrotron emission associated with the intra- 
cluster plasma medium. Recent lists of these objects can 
be found in, e.g., Feretti (1999), Giovannini (1999), Gio- 
vannini, Tordi & Feretti (1999), Owen, Morrison & Voges 
(1999), and Giovannini & Feretti (2000). The radio mor- 
phology is classified into radio halos and relics depend- 
ing on whether (almost unpolarized) radio emission is de- 
tected throughout the cluster, or (higher polarized) radio 
emission is detected in the cluster periphery. Radio halos 
and relics are predominantly found in rich compact clus- 
ters with Bautz-Morgan types II-III (Hanisch 1982), have 
high X-ray luminosities and temperatures, in most cases 
no cooling flow signature, and large cluster X-ray core 
radii. It is suggested (see, e.g., Harris et al. 1980, Burns 
et al. 1995, Feretti & Giovannini 1996, and the review of 
Sarazin 2001) that radio halos form during the merging of 
subclusters, accelerating (existing relativistic) particles in 
shocks formed in the intracluster gas, although additional 
processes seem to be necessary. Because of the compara- 
tively small sample sizes involved (see Sect.|^) we refer in 
the following only to the combined radio halo/relic sam- 



ple, neglecting the interesting differences between the two 
cluster types. 

Section^ describes the methods used to detect sub- 
structure in two-dimensional photon distributions and to 
assign statistical significances to the results. The REFLEX 
and BCS samples as well as the radio halo/relic and cool- 
ing flow samples are described in Sect.||. Further tests and 
illustrations of the morphological classifications are given 
in Sect.|I| The observed frequency distributions of clusters 
with substructure are discussed in Sect.||. The substruc- 
ture occurrence rates (SORs) show a comparatively strong 
sensitivity to the number of X-ray photons used for sub- 
structure analyses. A simple method to reduce this effect 
is described. 

To proceed further, two methods seem useful. The first 
possibility is the application of our substructure tests to 
simulated images where selection effects introduced by 
observation and data reduction are taken into account. 
In addition, the comparison with simulations is expected 
to yield statistical information about the mass spectrum 
of the merger masses to which our substructure tests are 
sensitive at various redshifts and X-ray fluxes. This work 
is postponed to a later paper. The second possibility is 
the relative comparison of subsamples extracted from the 
combined REFLEX+BCS catalogs in a manner that selec- 
tion effects partially cancel out. The subsequent sections 
of this paper are devoted to this kind of analysis. 

One important result obtained with the latter strat- 
egy is the discovery of a substructure-density relation, 
analogous to the morphology-density relation for galax- 
ies (Sect.||). Other interesting results are obtained from 
the relative comparison of substructure significance dis- 
tributions of halo/relic, cooling flow, and REFLEX+BCS 
clusters (Sect.^). Here we create reference subsamples 
drawn randomly from the REFLEX+BCS sample in such 
a way that the reference samples have basically the same 
'sensitivity' for substructure detection as the halo/relic 
and cooling flow samples. With this method we can com- 
pare for the first time the substructuring and elongation 
behaviour of specific cluster types with a representative 
large-area cluster sample. A general discussion of the re- 
sults is given in Sect.|[ 

All computations assume Friedmann-Lemaitre world 
models with zero pressure, the Hubble constant Hq = 
50 kms -1 Mpc -1 , the density parameter Slo = 1-0, and 
the cosmological constant = 0. 

2. Statistical tests for substructure 

2.1. General considerations 

In X-rays, quantitative analyses of substructure are based 
on detailed fits of elliptical models or wavelets to X-ray 
isophotes, studies of ellipticities, center-shifts, power ra- 
tios, etc. (see, e.g., the review in Buote 2001). However, 
representative SORs can only be obtained with large clus- 
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ter samples. The RASS-3 X-ray images of our cluster sam- 
ple which can be used for this purpose have on average 
359 with up to 6 829 X-ray photons (in the ROSAT hard 
energy band, see below) within an aperture of 1 Mpc ra- 
dius. Therefore, more robust tests for substructure are 
needed. Unfortunately, the relation between substructure 
as defined by robust tests and a physical quantification 
of substructure is less direct and a larger effort is needed 
for the physical interpretation of the results. The link be- 
tween SORs and theoretical merger rates can be obtained 
when the mass scales of the subclumps and the time scales 
needed for the merged cluster to reach dynamical equilib- 
rium are known. Many observational effects become ap- 
parent when statistical samples are analysed (see Sects. 0, 



5.1 



) which complicate the interpretation of SORs. After- 
proper correction, semi-analytic Press Schechter-like the- 
ories as presented in, e.g., Lacey & Cole (1993) can in 
principle be applied to understand the merger rates within 
a cosmological context. 

However, the relation between observed SORs and true 
substructure becomes secondary when SORs obtained for 
different subsamples drawn from the same parent distri- 
bution are compared. In this case the tests define sub- 
structure more operationally as that what they measure, 
and they can serve as a mere link between different clus- 
ter types as, for example, clusters located in high and low- 
density regions (Sect.|6|) or halo/relic clusters, cooling flow 
clusters (Sect.0). 

N-body simulations of merging clusters of galaxies 
favour three tests for the analysis of two-dimensional point 
distributions (Pinkney et al. 1996). The tests described in 
Sects. ^31 to 2.4 are sensitive to different types of substruc- 



ture and are thus ideally suited for the detection of a large 
variety of different merger events. We use the three tests 
and translate them to the case of two-dimensional X-ray 
images as extracted from the RASS-3 fields. Sect.[| com- 
pares substructures as defined in the present paper under 
realistic conditions with the results obtained by other re- 
search groups with different methods. 

2.2. X-ray survey data 

Before the substructure statistics are described the ba- 
sic properties of the X-ray material used for the analy- 
sis are summarized. This will explain some of the con- 
straints already imposed by the observational data. The 
1 378 RASS-3 fields used for substructure analyses cover 
the whole sky with an averaged spatial resolution as given 
by the half power radius of 96 arcsec at 1 keV. Substruc- 
ture on scales larger than approximately 2Mpc-z (e.g., 
200 kpc at z — 0.1) is thus resolved. Within this precision, 
the point spread function is constant over each RASS field 
and can be regarded as circular symmetric (G. Boese, pri- 
vate communication). Each field has a size of 6.4 x 6.4 deg 2 
and overlaps at least 0.23 deg with adjacent fields. The 
main advantage of RASS-3 compared to RASS-2 is that 



its less stringent constraints on the attitude solutions yield 
a larger number of accepted X-ray photon events (on aver- 
age about 5 percent, and is essential only in certain parts 
of the sky) resulting in a higher signal-to-noise without a 
significant increase of the measurement errors of the indi- 
vidual photons. The advantages of RASS-2 against RASS- 
1 are discussed in Voges et al. (1999). 

Each of the Af p h photons detected in the direction of a 
cluster of galaxies is characterized by the following quan- 
tities. The sky pixel coordinates, {xi, yi), given for the ith 
photon in units of 0.5 arcsec, and the photon energy, 
given in the form of a Pulse Height Amplitude (PHA) 
channel number, are contained in the photon event file of 
each RASS-3 field. The exposure time, ti, used to weight 
each photon, is determined with the exposure map, giving 
for each RASS-3 field the survey exposure times, corrected 
for vignetting and the effect of the shadowing of the sup- 
port structure of the PSPC window with a spatial binning 
of 45 arcsec. We do not differentiate between source and 
background photons because no distinction can be made 
on the photon-by-photon basis. 

Only X-ray photons selected from the ROSAT hard 
energy band (PHA channels 52 to 201 corresponding to 
about 0.5-2.0 keV) are used for all reductions. This reduces 
the soft X-ray background by a factor of 4, but still keeps 
60 to 100 percent (depending on the interstellar column 
density) of the cluster emission. Moreover, the majority 
of soft RASS-3 sources superposed onto the cluster im- 
ages are suppressed, so that the signal-to-noise for cluster 
detection is highest. 

Since comparatively large samples of clusters have to 
be analyzed the computer-intensive substructure tests had 
to be simplified (see below). The tests are performed 
within one aperture with the metric radius of 1.0 Mpc cen- 
tered on the X-ray intensity- weighted center of the cluster. 
The iterative determination of the cluster center termi- 
nates when a formal accuracy of 2 arcsec is reached (for 
more details see Bohringer et al. 2000, 2001). 

2.3. The (3 test 

The [3 statistic compares the surface number density, fa, 
measured in the ith radial segment with the density, \i i, 
measured in the diametrically opposite segment. The num- 
ber densities are obtained by weighting each photon with 
the corresponding exposure time. The final (3 value is the 
average of the ratios obtained over all independent N s ra- 
dial segments (typically 8): 

This test is a simplified version of a method developed by 
West, Oemler, & Dekel (1988). 

As noted in Pinkney et al. (1996) the (3 value (in its 
original definition) is sensitive to deviations from mirror 
symmetry, independent of the actual elongation of the tar- 
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get. The statistic becomes rather ineffective as the mass 
ratio of the sub-components approaches f : f . We expect 
similar properties of the presently implemented simplified 
version. Note that the application of radial segments in- 
stead of individual photon coordinates is used to make the 
test more robust and to save computing time although this 
might desensitize the test. The actual sensitivity of the (3 
parameter is tested with RASS-3 images in Sect.||. 

2.4- The Lee test 

In the Lee statistic the photons are projected onto lines 
with given inclination angles, (f> (12 directions in the 
present implementation). In general, for all <f> and for all 
partitions of the set of photons, the 'within-class' scatter, 
cl + or.) and the 'between-class' scatter, <tt, are deter- 
mined and their ratios, ot/^l + or), maximized. We use 
the test based on the technique of maximum likelihood in 
the form 

imax _ max^ {L((f>)} 



L = 



where 

L(4>) = max{ partitiolls } 



{L(4>)} 



0"X 



1 



(2) 



(3) 



, cr L + CTr 

gives the maximum likelihood obtained for the (-/V p h — 1) 
partitions. The partitions are obtained by dividing the 
total set of photon coordinates (T) projected onto a line 
with a given (j) into a right (R) and a left (L) subgroup 
with the corresponding scatters 

Ei _ - \2 - _ 2^j=i w j x j i .s 

Wj (Xj x a ) , x a — x > 

where a = T, R, L, the weighting factor for the ith photon 
is Wj = 1/tj, and K a is the actual number of projected 
points in the T, R, L samples. More details can be found 
in Fitchett et al. (1988). Fitchett & Webster (1987) suc- 
cessfully applied the method for substructure detection in 
the core of the Coma galaxy cluster. 

The likelihood value, L, is most sensitive if two sub- 
structure components are present, especially when they 
are very compact. It is not sensitive to any elongations 
and loses sensitivity if more than two subclumps are visi- 
ble. Compared to the other two tests the LEE statistic is 
thus the most conservative. 



2.5. The Fourier elongation test 

Deviations from circular symmetry might indicate merger 
events (e.g., Rocttigcr et al. 1997, see also Buote 2001), 
although relaxed clusters can show significant ellipticities 
too as illustrated in Sect.|]. However, recent simulations 
of Thomas et al. (1998) and others (see Sect.Q) clearly 
show a dependency of the frequency distribution of clus- 
ter major axial ratios on cosmology, where the rounder 
X-ray isophotes in the low-ilo models might result pri- 
marily from the scarcity of recent mergers (Evrard et al. 
1993). We thus regard elongation as a useful cosmological 



quantity although its relation to substructure appears to 
be more complicated. 

Following Rhce, van Haarlem & Katgert (1991) the az- 
imuthal number counts are approximated to first order by 
a constant density, modulated by a double sine. Under this 
assumption the normalized amplitude of this modulation, 



FEL 



'2(5 2 + C 2 ) 



(5) 



gives a measure of the elongation strength. The parame- 
ters S and C are defined by equation (7) in Pinkney et al. 
(1996). The position angle is then PA = 0.5 arctan(SyC). 
We translate (||) to the unbinned case, where each photon 
is weighted by its exposure time, 

JVph „ iVpfc 2 2 

with No = J2i The N-body simulations show that 
FEL is the most sensitive of the three tests used to check 
for substructure if elongation is to be considered as sub- 
structure. 

2.6. Statistical significances 

In the following the method is described which is used 
to compute the probabilities (statistical significances S) 
that the actual valu es of the substructure parameters de- 
scribed in Sects. 2J3 to 2A could be obtained just by chance 
from an X-ray image satisfying the null hypothesis of a 
circular symmetric, mirror-symmetric, and unimodal sur- 
face brightness distribution. The statistical significance, 
S, thus corresponds to the confidence probability, (1 — S), 
that the null hypothesis can be rejected. 

The significance values, Sp, Sfel, <Slee, are computed 
by comparing the /3, FEL, and L values obtained for 
the programme cluster with the corresponding values ob- 
tained with a large set of unstructured photon distribu- 
tions (replicants) derived from the same cluster. For the 
position-independent, circular symm etric average point 
spread function of RASS-3 (see Sect.|2.2|), these smooth 
and symmetric distributions can be obtained by azimuthal 
randomization (West et al. 1988), keeping the radial dis- 
tances of the original X-ray photons, but assigning random 
angular coordinates without the need to choose a specific 
(model) cluster profile. 

The significances are computed as follows. In the first 
step /3, FEL, and L are computed for the programme clus- 
ter. In the next step the photon distribution is azimuthally 
randomized, and for each randomized cluster, /3', FEL', 
and L' are computed again. The primes denote values 
of the statistics obtained with the randomized distribu- 
tion. As seen above the statistics are normalized in such 
a way that larger values correspond to larger substruc- 
tures. Therefore, the number of times the randomized dis- 
tribution gives a value larger than that obtained for the 
programme cluster provides an estimate of the probabil- 
ity, S, that the actual substructure value can be obtained 
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just by chance from a cluster fulfilling the null hypothesis. 
Notice that small values of S correspond to clusters with 
substructure and elongation. 

To be more specific, the substructure analyses are per- 
formed and statistical significances are computed for one 
aperture with the metric radius of 1.0 Mpc. This metric 
scale is transformed into an angular scale with the cluster 
redshift using the angular diameter distance. The choice 
of this radius guarantees that in the majority of cases the 
outer significance radii of the X-ray images, defined below, 
include the aperture used to evaluate substructure. For 
larger apertures, contributions from neighbouring clus- 
ters and chance superposition of background sources be- 
come more important. The significance radius is defined 
in Bohringer et al. (2000, 2001) as the point where the 
increase in the la flux error is larger than the increase of 
the cumulative source count rate after background sub- 
traction. For the (3 and FEL statistic, 2 000 replicants are 
used for each cluster, for the most time-consuming LEE 
statistic 400 replicants. 

3. Cluster samples 

The REFLEX sample (Bohringer et al. 2001) includes 
452 clusters (449 with redshifts), detected in the RASS- 
2 database (Voges et al. 1999) south of the Declina- 
tion +2.5 deg and down to the nominal X-ray flux limit 
of 3.0 x 10~ 12 ergs -1 cm~ 2 in the ROSAT energy band 
(0.1-2.4keV). REFLEX covers 4.24sr excluding the area 
±20 deg around the galactic plane and 0.0987 sr around 
the Magellanic Clouds and is basically restricted to red- 
shifts z < 0.3 (a few clusters reach z — 0.45). The Growth 
Curve Analysis method of Bohringer et al. (2000, 2001) 
is used to compute source positions, fluxes, angular ex- 
tents, etc. (source characterization). The catalog gives un- 
absorbed X-ray fluxes with statistical errors between 10 
and 20 percent. Several tests indicate a high overall com- 
pleteness of the sample (at least 90 percent) with an upper 
limit of 10 percent of the clusters with fluxes significantly 
contaminated by active galactic nuclei (AGN). 

The BCS sample (Ebeling et al. 1998) includes 201 
clusters (statistical sample), detected north of the Decli- 
nation 6 = Odeg, excluding the area ±20 deg around the 
galactic plane using the RASS-1 data. The formal flux 
limit of the resulting sample is 4.4 x 10 -12 erg s _1 cm -2 
(0.1-2.4keV). The Voronoi Tesselation and Percolation 
method (first version of the method described in Ebeling 
& Wiedenmann 1993) is used to obtain unabsorbed X-ray 
fluxes. The authors estimated a sample completeness of 
90 percent for the 201 BCS clusters with z < 0.3. Some 
artificial fluctuations in cluster number density might oc- 
cur because the clusters are not sampled homogeneously 
over the total survey area given in Ebeling et al. (1998). 
For the present investigation the effects can be neglected 
because only large-amplitude fluctuations (> 10 precent) 
are discussed. 



There are basically two motivations for using a com- 
bined sample of REFLEX and BCS clusters for our sub- 
structure analysis. (1) It is expected that the results of 
substructure analyses are improved when more X-ray pho- 
tons are available (see Sect.||). Combining the REFLEX 
and BCS catalogues to an all-sky sample and cutting it 
at the flux limit of the brighter BCS sample yields the 
largest number of X-ray images with high numbers of X- 
ray photons. (2) Due to the fact that at low redshifts the 
northern BCS sample might be dominated by the local su- 
percluster and its extension, and the southern REFLEX 
sample by the southern void (e.g., Schuecker et al. 2001), 
the merged REFLEX and BCS catalogue is expected to 
provide a more representative sample. 

To ensure an homogeneous handling of the data we 
exclude the Virgo cluster because of its large angular ex- 
tent. Furthermore, A689, Z3179, and RXCJ1212. 3-1817 
are excluded because of their low ROSAT exposure times 
(too few photons), A1678 and RXCJ0532.9-3701 which 
are located too close to the RASS-3 field edge, and three 
further REFLEX clusters where redshift information is 
not yet available. The combined sample referred to as 
REFLEX+BCS contains 470 clusters (excluding the BCS 
clusters in common with REFLEX) with a formal flux 
limit of 4.4 x 10~ 12 erg s" 1 cm" 2 . Although BCS and RE- 
FLEX fluxes are computed with different algorithms on 
different RASS versions no problems are expected because 
each version has its own calibrations. 

A sample of 53 halo and relic clusters (including six 
uncertain cases) compiled by Feretti et al. (in prepara- 
tion) is used to study the statistical properties of their 
X-ray morphology. In the following we use the combined 
halo/relic sample as mentioned in Sect. |l|. Deep pointings 
obtained mainly with ROSAT suggest that eventually all 
of these clusters show distorted X-ray surface brightness 
distributions. Presently no information about the com- 
pleteness of this sample is available. In the present inves- 
tigation RASS-3 data are used for their analysis so that 
the substructure tests can be directly compared to the sta- 
tistically representative REFLEX+BCS reference sample. 

A sample of 22 clusters with large cooling flow sig- 
natures is selected from the list presented in Peres et 
al. (1998). Their total sample of 55 clusters with fluxes 
above 1.7x 10 -11 ergs -1 cm -2 in the 2-10 keV energy band 
is selected from observations with the Einstein and EX- 
OSAT observatories, and the HEAO-1 and Ariel V satel- 
lites by Edge et al. (1990) and is found to be 'satisfactorily 
complete'. The selected cooling flow clusters have mass 
deposition rates larger than 100A/©yr -1 as determined 
by a surface brightness deprojection technique using RO- 
SAT PSPC and HRI pointed observations. Notice that for 
the present investigation only the cooling flow signature 
of a steep increase of the central X-ray surface brightness 
matters, and that a possibly new interpretation of cooling 
flows due to the non-detection of cooling gas with tem- 
peratures below about 3keV in XMM spectra (Peterson 
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Fig. 1. Comparison of iso-surface brightness contours of deep ROSAT PSPC pointings (first and third columns) and corre- 
sponding RASS-3 images (second and fourth columns) of regular clusters. The pointings are superposed on Digital Sky Survey 
optical images. Source characterizations are given in Tab. El (clusters 1-10). 



et al. 2001) is secondary. The majority of the clusters in 
the sample have z < 0.1. 

4. Verification of the method 

The methods used for substructure detection were ana- 
lyzed and tested in detail with numerical N-body simu- 
lations (Pinkney et al. 1996). Here we want to illustrate 
the types of substructures to which the three substructure 



tests are sensitive under realistic conditions comparing 
deep ROSAT pointed observations and the correspond- 
ing RASS-3 images of clusters with known morphology 
and substructure behaviour (see Sect. 4.1). In addition, 
we want to study the relation between the substructure 
tests and the Jones & Forman classification scheme to al- 
low a reference to a similar systematic survey project (see 
Sect. 4.2). It will be seen that the three substructure statis- 
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Fig. 2. Comparison of iso-surface brightness contours of deep ROSAT PSPC pointings (first and third columns) and corre- 
sponding RASS-3 images (second and fourth columns) of substructured clusters. The pointings are superposed on Digital Sky 
Survey optical images. Source characterizations are given in Tab. | (clusters 11-20). 



tics provide a quantitative morphological classification of 
RASS-3 images which has a close relation to substructure 
as defined and quantified with other methods and data. 

4-1- Comparison with deep ROSAT PSPC pointings 

Figures and | show a collection of 10 regular and 10 
substructured 'prototype' clusters as obtained with deep 



ROSAT PSPC pointings (1st and 3rd columns) and the 
corresponding RASS-3 images (2nd and 4th columns) used 
for our substructure analyses. Literature information and 
the significances obtained with our substructure tests us- 
ing RASS-3 data within a metric aperture of 1.0 Mpc 
are summarized in Tabs.[j] and ||. From the literature we 
summarize the classifications obtained with the power ra- 
tio technique from ROSAT PSPC pointed observation by 
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Table 1. Characterization of the regular galaxy clusters shown in Fig. |T|. The mass deposition rates are in units of 
solar mass per year as obtained from ROSAT PSPC/HRI pointings, .Slee, Sp, 5fel are the statistical significances 
obtained with the RASS-3 images within circles of radius R in arcsec corresponding to a metric radius of 1.0 Mpc, 
position angles, PA in degrees as obtained with FEL and surface brightness-weighted moments. 



Cluster Characterization Reference 



1 A478 

2 A1651 

3 A1795 

4 A2029 

5 A2204 

6 A2244 

7 A2597 

8 A4059 

9 2A0335 
10 Hydra A 



smooth 

no centroid shift 

M = 520 - 616 

R = 8': S'lee = 0.28, Sg = 

mostly smooth+ regular 

M = 138 

R = 8': S'lee = 0.56, Sg = 
smooth + regular 
M = 381/488 

R= 10': S'lee = 0.15, S = 

regular + smooth 

M = 554/556 

R = 9': Slee = 0.87, Sg = 

smooth 

M = 843/852 

7? = 5': Slee = 0.59, Sg = 

regular 

M = 244 

R = 7': Slee = 0.13, Sg = 
very symm. + smooth 
M = 276/271 

R = 8': Slee = 0.16, Sg = 
smooth, single-comp. 
M = 130 

7?= 13': Slee = 0.96, Sg = 
M = 242/325 

7?= 18': Slee = 0.65, Sg = 

regular 

M = 298/264 

7?= 12': Slee = 0.35, Sg = 



0.76, Seel = 0.00, PA = 58/57 

0.47, Seel = 0.53, PA = 100/105 

= 0.64, Seel = 0.03, PA = 20/12 

0.07, Seel = 0.81, PA = 169/176 

0.45, Seel = 0.17, PA = 77/92 

0.10, Seel = 0.02, PA = 6/174 

0.49, Seel = 0.02, PA = 93/112 

= 0.30, Seel = 0.01, PA = 127/130 
= 0.06, Seel = 0.52, PA = 132/146 

= 0.37, Seel = 0.66, PA = 140/126 



Buote & Tsai (1996) 
Mohr et al. (1995) 
Peres et al. (1998) 

Buote & Tsai (1996) 
Peres et al. (1998) 

Buote & Tsai (1996) 
Peres et al. (1998) 

Buote & Tsai (1996) 
Peres et al. (1998) 

Buote & Tsai (1996) 
Peres et al. (1998) 



Buote & 
Peres et 

Buote & 
Peres et 

Buote & 
Peres et 



Tsai (1996) 
al. (1998) 

Tsai (1996) 
al. (1998) 

Tsai (1996) 
al. (1998) 



Peres et al. (1998) 



Buote & 
Peres et 



Tsai (1996) 
al. (1998) 



Buote & Tsai (1996), cooling flow mass deposition rates 
and cooling times as derived by Peres et al. (1998) us- 
ing ROSAT PSPC and HRI data, centroid variances de- 
duced by Mohr et al. (1995), and visual classifications 
of Jones & Forman (1999), both obtained with Einstein 
IPC data, ellipticities and center-of-mass shifts as de- 
duced from Abell+APM galaxy and ROSAT X-ray data 
by Kolokotronis et al. (2000), and results obtained with 
the adaptive-kernel technique by Kricsslcr & Beers (1997) 
using optical galaxies in Dressler's clusters. We also add 
the presence of cold fronts as given in Markevitch et al. 
(2000) and Vikhlinin et al. (2001). 

The clusters shown in Fig.|l| are nearby z < 0.1 clus- 
ters, all with strong cooling flow signatures and classified 
as single, regular, smooth, and/or no centroid shift by the 
authors given in Tabs.|l] and [| The S'lee significances 
as obtained from the RASS-3 images range between 0.13 
and 0.96, the Sp values between 0.06 and 0.76, and the 
Sfel values between 0.0 and 0.81. Assuming a 99 percent 
confidence level none of the LEE and f3 tests indicate sig- 



nificant substructure for these clusters. According to our 
substructure statistics, clusters like A478 and A4059 show 
a significant ellipticity, including A1795, A2244, A2597 on 
the 97 percent confidence level (in total 5 out of the 10 
clusters). If we assume that the central 1.0 Mpc of these 
clusters, where the tests are performed, are dynamically 
relaxed and unstructured then ellipticity and thus FEL 
would turn out to be a less useful indicator for substruc- 
ture whereas j3 and LEE would do. 

The non-regular clusters shown in Fig.^ have negligi- 
ble cooling flows, some have cold fronts. The clusters are 
classified as substructured by different authors using vari- 
ous methods and data (X-ray, optical). For A754, A1367, 
and A3266 the f3, LEE, and FEL tests all give zero proba- 
bility for the null hypothesis, thus strongly suggesting sub- 
structured RASS-3 surface brightness distributions. For 
A119, A2256, and A3667 either the (3 and/or the LEE 
test suggest the presence of substructure. On the 98 per- 
cent confidence level this is also true for A400, A2657, 
A3562, and A3921. Therefore, basically all 10 prototype 
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Table 2. Characterization of the substructured galaxy clusters shown in Fig.g. The mass deposition rates are in units 
of solar mass per year as obtained from ROSAT PSPC/HRI pointings, Slee, Sp, Sfel are the statistical significances 
obtained with the RASS-3 images within circles of radius R in arcsec corresponding to a metric radius of 1.0 Mpc, 
position angles PA in degrees as obtained with FEL and surface brightness-weighted moments. 



Cluster 


Characterization 


Reference 












11 


A119 


emission tail 










Buote & Tsai (1996) 






opt. substructured 










Kriessler & Beers (1997) 






centroid shift 










Mohr et al. (1995) 






M = 










Peres et al. (1998) 






R= 14': Slee = 0.84, 


Sg 

P 


= 0.00, Sfel 


= 0.13, PA 


= 17/27 




12 


A400 


irregular 










Buote & Tsai (1996) 






centroid shift 










Mohr et al. (1995) 






R = 25': Slee = 0.02, 


Sg 


= 0.34, Sfel 


= 0.11, PA 


= 176/155 




13 


A754 


clearly off center 










Buote & Tsai (1996) 






opt. substructured 










Kriessler & Beers (1997) 






centroid shift 










Mohr et al. (1995) 






M = 0/2 










Peres et al. (1998) 






Radio halo + relic 










Cohen et al. (2001) 






Cold Front 










Markevitch et al. 2001 






R = 12': Slee = 0.00, 


Sg 


= 0.00, Sfel 


= 0.00, PA 


= 101/101 




14 


A1367 


centroid shift 










Mohr et al. (1995) 






M = 8/0 










Peres et al. (1998) 






R = 28': Slee = 0.00, 


Sb 


= 0.00, Sfel 


= 0.00, PA 


= 146/145 




15 


A2256 


unrelaxed 










Buote & Tsai (1996) 






elliptical 










Jones & Forman (1999) 






substructured core 










Mohr et al. (1995) 






M = 0/0 










Peres et al. (1998) 






Radio halo + relic 










Feretti (1999) 






R = 11': Slee = 0.00, 




= 0.03, Sfel 


= 0.00, PA 


= 128/128 




16 


A2657 


bi-modal 










Buote & Tsai (1996) 






opt. substructured 










Kriessler & Beers (1997) 






R= 15': Slee = 0.02, 


Sg 


= 0.23, Sfel 


= 0.00, PA 


= 77/85 




17 


A3266 


substructured core 










Mohr et al. (1995) 






strong substructure 










Kolokotronis et al. (2000) 






M = 0/3 










Peres et al. (1998) 






R = 11': Slee = 0.00, 


Sp 


= 0.00, Sfel 


= 0.00, PA 


= 65/63 




18 


A3562 


distorted/not ellip. 










Buote & Tsai (1996) 






M = 37 










Peres et al. (1998) 






R= 12': Slee = 0.68, 


Sp 


= 0.02, Sfel 


= 0.60, PA 


= 94/61 




19 


A3667 


elongated + distort. 










Buote & Tsai (1996) 






apparent bi-modal 










Kolokotronis et al. (2000) 






M = 










Peres et al. (1998) 






Cold Front 










Vikhlinin et al. (2001) 






2 radio relics 










Feretti (1999) 






R = 11': Slee = 0.00, 


Sp 


= 0.80, Sfel 


= 0.00, PA 


= 133/113 




20 


A3921 


irregular/merger 










Buote & Tsai (1996) 






apparent bi-modal 










Kolokotronis et al. (2000) 






R=7': Slee = 0.02, Sp = 


= 0.88, Sfel = 


= 0.00, PA = 


= 99/102 





clusters shown here are recovered as substructured on the 
98 percent confidence level with the (3 and/or LEE test on 
RASS-3. 

We found a good correspondence between the position 
angles, PA, as seen in the deep pointings and as obtained 
from the RASS-3 images (given in Tabs.|l] and |). The 
position angles are obtained with both the FEL method 
and with surface brightness- weighted moments. In basi- 
cally all cases with elongation significances Sfel < 0.01, 



position angles obtained with the FEL method and those 
determined with surface brightness-weighted moments are 
found to be very similar within a few degrees. 

Moreover, 7 of the 10 clusters (A754, A1367, A2256, 
A2657, A3266, A3667, A3921) have a significant elliptic- 
ity. This suggests a correlation between substructure and 
the detection of ellipticity. However, as the analysis of the 
cooling flow clusters already showed, ellipticity as a sole 
criterion for substructure might not be sufficient for sub- 
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Table 3. Selected HIFLUGCS clusters and their statistical significances for the LEE, /3, and FEL statistics for regular 
(Reg), possibly substructured (Poss), and substructured (Sub) clusters. The significances are obtained from RASS-3 
images, the classification of the clusters is taken from the literature. Mean significances (Meanl) and their formal 
standard deviation (Stdl) as well as the mean significances (Mean2) and their standard deviation (Std2) obtained for 
the minimum values of (Slee, Sp) are given in the last rows (eq. ^). Notice that A85 has a substructure just outside 
the 1 Mpc aperature and is thus regarded as possibly substructured. 





Ree- 








Poss 








Sub 






Name 


5*LEE 


Up 


*JFEL 


Name 


Slee 


S i 


q 

^FEL 


Name 


•Slee 


S 4 

p 


^FEL 


A 478 


n 98 


n 7fi 

U. I u 


n nn 

u.uu 


AS^ 


n 48 


0.21 


n on 


A119 


n 84 


n nn 

u.uu 


n 1 3 

U. -LO 


A1651 


U. UU 


n 47 


n ^ 


A 1 11 


n to 


U. -L 1 


n 1 1 

u. ±o 


A400 


n D9 


n id 


n 1 1 


A 1 70^ 


n i <\ 


n is4 


n m 

U. uo 


A401 


n 07 

u.u i 


0.22 


VJ.XJO 


A496 


n m 


0.14 


n nn 

u.uu 


A909Q 


0.87 


0.07 


0.81 


A 9063 


0.70 


0.01 


0.77 


A754 


0.00 


0.00 


0.00 


A2052 


0.80 


0.02 


0.06 


A3158 


0.01 


0.13 


0.00 


A1367 


0.02 


0.00 


0.00 


A2204 


0.59 


0.45 


0.17 


A 3558 


0.07 


0.06 


0.00 


A2142 


0.04 


0.00 


0.05 


A2244 


0.13 


0.10 


0.02 










A2147 


0.00 


0.00 


0.00 


A 2589 


0.01 


0.39 


0.00 










A 2255 


0.00 


0.00 


0.00 


A2597 


0.16 


0.49 


0.02 










A2256 


0.00 


0.03 


0.00 


A3112 


0.39 


0.90 


0.02 










A2634 


0.02 


0.01 


0.02 


A3581 


0.13 


0.96 


0.13 










A2657 


0.02 


0.23 


0.00 


A4038 


0.00 


0.95 


0.00 










A3266 


0.00 


0.00 


0.00 


A4059 


0.96 


0.30 


0.01 










A3376 


0.00 


0.00 


0.00 


MKW3s 


0.24 


0.44 


0.21 










A3391 


0.00 


0.46 


0.00 


Hydra-A 


0.35 


0.37 


0.66 










A3395 


0.00 


0.00 


0.00 


S1101 


0.73 


0.11 


0.62 










A3562 


0.68 


0.02 


0.60 


2A0335 


0.65 


0.06 


0.52 










A3571 


0.00 


0.92 


0.00 


















A3667 


0.00 


0.80 


0.00 


















A3921 


0.02 


0.88 


0.00 


















RXJ2344 


0.04 


0.02 


0.06 


Meanl 


0.412 


0.440 


0.224 




0.320 


0.123 


0.155 




0.086 


0.193 


0.049 


Stdl 


0.075 


0.077 


0.068 




0.124 


0.034 


0.125 




0.052 


0.071 


0.030 


Mean2 


0.194 








0.078 








0.006 






Std2 


0.038 








0.031 








0.002 







structure detection. Nevertheless a final decision whether 
ellipticity is regarded as substructure is not necessary for 
the present investigation. 

In addition to the 20 prototype clusters discussed 
above, Tab.|| summarizes the significances obtained for 
43 clusters included in the Highest X-ray FLUx Galaxy 
Cluster Sample (HIFLUGCS, see the sample description 
in Reiprich & Bohringer 2001) which have ROSAT PSPC 
pointed observations and reliable literature information 
of the kind given in Tabs.[l] and || about the presence or 
absence of substructure (see the references given at the be- 
ginning of this section). The mean significances (Meanl) 
and their formal ler deviations (Stdl) given in Tab. | show 
a clear trend towards smaller significances for the clusters 
expected to be substructured. This trend is even stronger 
when mean and standard deviation are computed for the 
minimum values 

S = min {Slee, Sp } , (7) 
as given by the Mean2 and Std2 values in Tab.|^. The last 
statistic (0) assumes that either a low Slee or Sp value 
already indicates substructure which seems to be reason- 
able in the light of the discussion of individual clusters 
given above. 



4-2. Comparison with Jones & Forman classifications 

Jones & Forman (1999) used the iso-intensity contour 
plots of the Einstein IPC X-ray emission of 208 targeted 
and serendipitously found clusters with z < 0.15 to clas- 
sify their morphology into the following categories. Single: 
no substructure or departures from symmetry. Double: two 
subclusters of comparable size and luminosity. Primary 
with small secondary: main subcluster at least two times 
brighter than secondary. Complex: more than two subclus- 
ters. Elliptical: elliptical X-ray surface brightness contours. 
Off center, peak emission not in center defined by lower 
surface brightness emission. Galaxy: emission dominated 
by a galaxy. The Einstein and RASS-3 X-ray images have 
similar angular resolutions (about 1 arcmin). Moreover, 
the Einstein classifications are restricted to clusters with 
redshifts z < 0.15. Therefore, one should expect some cor- 
relation between the Jones & Forman and our classifica- 
tions. 

For the comparison of the two classification schemes we 
follow first Jones & Forman and decide that only the Sin- 
gle and Galaxy types are regarded as regular and almost 
relaxed structures. The remaining classes (Double, Pri- 
mary with small secondary, Complex, Elliptical, Off cen- 
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Fig. 3. Average significances and their standard deviation as a 
function of the number of RASS-3 photons of the 116 Einstein 
clusters in common with the BCS and REFLEX cluster sam- 
ples. The upper curves in each of the three panels is obtained 
for the Einstein clusters classified as single by Jones & Forman 
(1999), the lower curves for the Einstein clusters classified as 
double, primary with small secondary, complex, elliptical, or 
offset center. 



ter) are regarded as clusters with significant substructure. 
At this stage we do not doubt the classification of Jones 
& Forman because we only want to illustrate the relation 
between two classification schemes. 

In total 116 classified Einstein clusters are found in 
common with the REFLEX and BCS sample and are re- 
garded as a representative test sample. Let us first show 
the relation of the three statistical significances with the 
classifications of Jones & Forman without introducing an 
ad hoc threshold for the confidence limit. Fig.|| compares 
the average significances obtained with each of the three 
statistics for the clusters classified by Jones & Forman as 
regular (upper continuous lines including their formal la 
Poisson errors) and substructured (lower dashed curves). 
The average significances are computed for subsamples 
as a function of the minimum number of X-ray photons, 
N m ia, m the RASS-3 images. 

For all N m i n values a clear trend is seen that clusters 
classified by Jones & Forman as regular have on average 
larger significances compared to the clusters classified by 
Jones & Forman as substructured. The largest differences 
between the mean S values of regular and substructured 
clusters are seen for the elongation statistic. This indicates 
that elongation or ellipticity is comparatively easy to clas- 
sify both by eye and machine, especially for large N^m- 



0.01 



0.001 




Fig. 4. Distribution of LEE significances as a function of the 
number of X-ray photons for REFELX+BCS clusters (points), 
for the substructured cluster A2256 (lower continuous line), 
and for the cooling flow cluster A2029 (upper continuous line). 
The lines represent significances obtained with different num- 
bers of X-ray photons selected randomly from the X-ray im- 
ages. Zero significance values are set to ~ 5 x 10 -4 . 



Similar results su ppo rting the robustness of this criterion 



are given in Sect. 5T 



In the next step we want to illustrate the conditions 
under which our substructure tests give a zero-order quan- 
titative implementation of the Jones & Forman morpho- 
logical classification scheme although the X-ray data used 
for the two sets of classifications are different. This also 
includes the introduction of a fixed confidence level for 
substructure classification. Consistent with the results ob- 



tained in Sect. 4.1 we now assume that clusters which are 
classified by Jones & Forman as elliptical are in fact un- 
structured. In order to get higher signal-to-noise surface 
brightness distributions and thus more reliable substruc- 
ture detections, we restrict the test sample to those clus- 
ters which have at least 200 X-ray photons on the corre- 
sponding RASS-3 images. This leaves 73 clusters for the 
final comparison. 

From this test sample, 59 are classified as regular and 
14 as substructured by Jones & Forman. We compare 
these numbers with the results obtained with our sub- 
structure tests assuming that clusters fulfilling the crite- 
rion min{5LEE, Sp} < 0.05 are substructured. In this case 
the substructure tests found 46 of the 59 clusters as regu- 
lar and 11 of the 14 clusters as substructured. This gives a 
fraction of 22 percent of the clusters with different classi- 
fications. The rate is slightly improved to 21 percent when 
a significance threshold of 0.01 is used. 

Notice that the values of the SORs depend on the de- 
tails of the sample properties as well as on the actual sub- 
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structure detection method so that an exact coincidence of 
SORs obtained with different samples cannot be expected. 
We conclude that the (3 and LEE statistic are found to 
provide useful tools for substructure analysis of RASS-3 
X-ray images and can under reasonable conditions approx- 
imately verify the Jones & Forman classification scheme. 

5. Occurrence rates of clusters with substructure 

The observed or apparent frequency of clusters with sig- 
nificant X-ray substructure is the basic quantity of our 
systematic study which can be determined with high accu- 
racy. The frequency distributions give useful information 
about the role of substructure and merger events in the 
nearby Universe, and provide several hints to important 
selection effects which might bias substructure statistics 
in general. 

For given spatial locations and size of main and sub- 
clusters, and for given angular resolution of the X-ray in- 
strument and aperture size, the detectability of substruc- 
ture depends on the number of X-ray photons, N p h, used 
to trace the cluster surface brightness distribution, and 
on the redshift of the cluster. Fig.^ illustrates the AT p h 
effect, plotting the statistical significances for the RE- 
FLEX+BCS clusters (dots) as a function of N ph . The 
superposed continuous lines give sequences of statistical 
significances obtained for the Abell clusters A2256 and 
A2029 for different numbers of photons selected randomly 
from the original RASS-3 photon distribution. A2256 is a 
well-known example of a cluster merger exhibiting a ra- 
dio halo and relic. Two subclumps are well-separated in 
the deeper ROSAT PSPC pointing (Briel et al. 1991, sup- 
ported by recent CHANDRA results given in Sun et al. 
2001), one maximum centered near the central cD galaxy. 
The cluster is located close to the North Ecliptic Pole 
so the correspondingly long ROSAT exposure time gives 
the comparatively large number of 3 000 X-ray photons 
within a circle of 1 Mpc in RASS-3 (hard band, includ- 
ing background photons). A2029 is a cooling flow cluster 
with a mass deposition rate > 550 solar masses per year 
and a cooling time < 3Gyr (Peres et al. 1998). Buote & 
Tsai (1996) classified this cluster as regular and smooth 
and used it as a prototype relaxed single cluster. Within 
a circle of 1 Mpc about 1 100 X-ray photons are found in 
RASS-3. 

For other clusters, and the (3 and FEL statistic, similar 
sample paths as shown by the continuous lines in Fig. ^ are 
found. Therefore, if clusters are located in the upper left 
part of the S-Nph diagram with say S > 0.01 correspond- 
ing to the 99 percent confidence level to reject the null hy- 
pothesis, this does not necessarily mean that they do not 
have any substructure. It only shows that under the actual 
detection conditions no statistically significant statement 
about the tested substructure property can be made. If 
more photons are collected the sample paths move to the 
right of the S-Np^ diagram. But whether the path moves 




Fig. 6. Significance curves for REFLEX+BCS clusters with 
redshifts z < 0.08 (dashed lines) compared to the template 
samples (continuous lines) covering the same redshift range. 
The input substructure fractions of the control sample used to 
fit the observations are shown in the upper left of each panel. 
The central continuous line represents the best fit of the sub- 
structure fraction to the observed curves. All curves are com- 
puted for the significance threshold S = 0.01. 

up or down would depend on the actual merger situation, 
the size of the subclusters, and on the cluster's redshift. 

The observed occurrence rates, /, of the 
REFLEX+BCS clusters obtained with the three sub- 
structure statistics are plotted in Fig.|| (left). For the 
morphological classifications the significance thresholds 
S < 0.02, 0.01, and 0.005 are applied. Independent 
of the corresponding confidence limits the observed 
fraction of substructured clusters increases with the 
minimum number of X-ray photons, N m { n . The compar- 
ison of the corresponding diagrams for subsamples of 
REFLEX+BCS clusters including an upper redshift limit 
(right panels of Fig.||) shows that redshift-dependent 
effects become important for minimum photon numbers 
A+mn < 500. They increase when in addition to an upper 
redshift limit also a lower limit is introduced. However, 
the fraction of high-z clusters is small so that their effect 
on the cumulative distributions is not large. Within the 
formal la Poisson errors, the observed fractions (complete 
sample) range between the lower limits of about 10 + 2 
percent and the plateau values of 60 ± 17 percent for the 
[3 and LEE statistic, and 78 ± 22 percent for the FEL 
statistic. 

The amplitude and shape of the curves depend on the 
'true' SOR and the X-ray flux and redshift distributions 
of the sample clusters which are in the end determined by 
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Fig. 5. Left panel: Fraction, /, of clusters as a function of the minimum number of X-ray photons, N m i n , with significances 
less or equal to the S = 0.005 (lower curves with formal Poisson error bars), 0.01 (middle curves) and 0.02 (upper curves) level. 
In the upper panel the measurements for the S = 0.01 and 0.02 levels degenerate above iV m m = 1 400 whereas for the lower 
panel all curves degenerate above -/V m j n = 800. The dashed horizontal lines mark the plateau values obtained for clusters with 
large numbers of X-ray photons. Right panel: Same curves as left panel for S — 0.01 and the redshift limits z — 0.04 (upper 
curves in each panel, < 88 clusters), z = 0.08 (< 239 clusters) and 0.45 (< 470 REFLEX+BCS clusters). 



the structure formation process, on the X-ray telescope 
and detector angular resolution, on the actual confidence 
level used for the classification, etc. 

5.1. Towards unbiased substructure occurrence rates 

Less biased SORs of the REFLEX+BCS clusters can be 
determined with a kind of template matching procedure 
where observed and model curves similar to those plot- 
ted in Figs.^J are compared and their differences mini- 
mized. The model curves are computed with a set of high 
signal-to-noise template clusters. This template set has by 
construction a known substructure occurrence rate, /inp, 
which can be changed by replacing substructured with 
regular templates or vice versa. Partially dependent repli- 
cants of these templates with different N-p^ are generated 
by diluting their X-ray images randomly in the same way 
as shown by the continuous sample paths in Fig.||. 

For a proper comparison the enlarged set of template 
images generated by dilution must have the same prob- 
ability density, P(N ph ), as the REFLEX+BCS clusters. 
Precise comparisons would also equalize the redshift distri- 
butions of observed and template clusters. This requires, 
however, a large set of high signal-to-noise X-ray images 
with at least 1000 X-ray photons at various redshifts. Such 
a template set is not available yet so that we simply re- 
strict the redshift range to be the same for both observed 
and template clusters. 



The matching procedure itself consists for each sub- 
structure criterion of the following steps. In the first 
step the fractions f(N m - ln ) are computed with the 
REFLEX+BCS sample (see Fig.|, left) and for a fixed 
confidence limit. In the second step an input substructure 
occurrence rate, /inp, of the template sample is chosen. 
The substructure test is computed for the diluted X-ray 
images. The images are randomly selected in a way that 
observed and template sets have the same P(N p \ 1 ) distri- 
butions. In the last step the curve, f(N m i n ), is determined 
and compared with the observed curve. Steps two and 
three are iterated by changing /inp of the templates until 
a useful fit of the observed SOR curve by the template 
model curve is achieved. 

Note that the application of a template sample based 
on RASS-3 data guarantees that many technical selection 
effects which might affect the REFLEX+BCS measure- 
ments are taken automatically into account. The quality 
of the results depends on the size and representativity of 
the template sample. 

The results shown in Fig.^| are computed for a tem- 
plate sample of 37 clusters (iV p h > 1000) with redshifts 
z < 0.08 and different input substructure fractions. No- 
tice that the classification of templates into substructured 
and regular clusters is done by using literature informa- 
tion of the kind given in Tabs.|j] and ||, and by visual 
inspection. As an example, curves are shown for a signif- 
icance threshold of S = 0.01. The results obtained with 
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the 239 RELFEX+BCS clusters are limited to the same 
maximum redshift as the control sample. Each template 
cluster is diluted 20 times so that an effective number of 
740 partially dependent template clusters is used. 

A simultaneous fit of the observed and model curves 
turns out to be complicated because the shapes of the 
curves differ systematically when specific iV m i n ranges 
are considered, indicating the presence of some remain- 
ing systematic differences between observed and template 
samples. One could, for example, expect that the tem- 
plate sample might not be representative of the complete 
REFLEX+BCS sample. Best (eye-ball) fits (continuous 
lines) give SORs ranging from 46 percent (FEL) to 59 
percent (/3). 

The reference (dotted) lines give some information 
about the sensitivity of the results to changes of /inp- 
They are computed for /inp = 0.39 and 0.75. The flat- 
ter FEL curves suggest that elongation is less affected by 
TVph. This result sup port s our evaluation of FEL for indi- 
vidual clusters (Sect. |4~2| ) and is quite important for our 
planned studies of alignment effects of position angles of 
major cluster axes. 

The final estimate of the 'true' SOR is estimated by 
the formal mean and standard deviation of the three input 
rates and gives 

SOR = 52 ± 7 percent . (8) 
Note that this result is independent whether elongation is 
regarded as substructure because FEL is only used as a 
link between template and REFLEX+BCS clusters. 

The nature of the merger events involved to yield the 
observed SOR is not clear because the applied substruc- 
ture tests do not give mass estimates of the individual 
subclumps. Therefore, numerical simulations are in prepa- 
ration to decide statistically which merger type (major 
merger or accretion) contributes at a given redshift and 
flux to the observed SOR. 

Another way to proceed further is the relative com- 
parison of subsamples derived from the REFLEX+BCS 
sample so that selection effects introduced by the actual 
observing conditions partially cancel out. Results obtained 
with this approach are described in the following sections. 

6. Substructure density relation 

Clusters in dense (supercluster) environments are ex- 
pected to have a higher probability to interact with neigh- 
bouring clusters or filamentary structures connecting the 
cluster centers. If this hypothesis is correct, larger frac- 
tions of clusters with distorted X-ray surface brightness 
distributions and thus with subclusters are expected in 
dense environments. 

As an illustrative example, Fig. [?] shows the spatial dis- 
tribution of the clusters with minlS^, Slee} < 0.1 (filled 
symbols) and > 0.1 (open symbols) in a region located at 
the Shapley concentration. It is seen that clusters with a 
lower significance to have a regular RASS-3 image (filled 
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-40° < 5 < -10° 

Fig. 7. Spatial distribution of clusters with min{S^, 5lee} < 
0.1 (filled symbols) and > 0.1 (open symbols). The radial coor- 
dinate is in kms -1 . For substructure detection a metric radius 
of 1 Mpc is used. 



symbols) appear to be preferentially more concentrated 
towards the core of the Shapley supercluster. If the sig- 
nificance threshold of 0.1 is lowered to 0.01 basically all 
clusters expected to show substructure are located in the 
Shapley concentration (only one cluster lies outside the su- 
percluster). The effect thus seems to be not very sensitive 
to the actual value of the significance threshold applied. 
Similar but less extreme trends are found for other super- 
cluster regions. 

To quantify this effect and to test its statistical sig- 
nificance, mean significances of (3 and LEE (for complete- 
ness we also use FEL) are computed for different local 
cluster number densities. For volume-limited samples, the 
computation of local cluster number densities would be 
straightforward. The REFLEX and BCS cluster samples 
used for the present investigation are, however, X-ray flux- 
limited and, although the samples are the largest yet avail- 
able, volume-limited subsamples derived from them have 
sample sizes smaller than about 100 (see, e.g., Tab. 1 in 
Schuecker et al. 2001). If we would further subdivide these 
subsamples with respect to different morphological classes 
the results would immediately suffer from small-number 
statistics. Therefore, we decided to work with flux-limited 
samples although the definition of rcdshift-independent 
measures of local cluster number density becomes less 
well-defined. 

As a measure of the local cluster number density 
around each cluster we use rf -3 , where d is the mean of its 
5 nearest neighbour distances. In order to correct for the 
redshift dependence of the cluster number density in flux- 
limited cluster samples, we normalize this number density 
by the average density obtained with the same density es- 
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Fig. 8. Left: Average significances and their standard deviation for the three statistical tests as a function of the local cluster 
number density, normalized to the average density for an aperture radius of 1 Mpc. Right: Same as left for an aperture radius 
of 3 Mpc. 



timator using all clusters in a redshift shell (z — Az, z+Az) 
centered on the cluster's z value (see comments given in 
Sect.|J). The normalization has the additional effect of 
compensating also for edge effects which are known to 
distort next neighbour statistics (e.g., Cressie 1993). 

The flux limit of REFLEX+BCS is quite bright and 
the segregation effect quite large (see below) so that it is 
not necessary to correct for local variations of the RE- 
FLEX and BCS survey sensitivity (e.g., RASS exposure 
convolved with galactic extinction). Note that cluster X- 
ray fluxes are, however, corrected for galactic extinction. 
Furthermore, no corrections for chance contamination of 
cluster images by non-cluster sources are applied because 
results obtained with different cluster subsamples drawn 
from the same parent distribution are compared. 

For the analysis of the REFLEX+BCS sample the av- 
erage cluster number densities, n, are computed within 
redshift shells with Az = 0.01 centered at the clus- 
ter redshifts. The local density contrasts are compared 
with the centers of optical superclusters given in Zucca 
et al. (1993). Among others, the Perseus-Pisces, Corona- 
Borcalis, Hercules, Shapley, and Horologium superclusters 
are clearly detected as regions with n/h > 1 for the ma- 
jority of the supercluster members. 

Figure || (left) shows the average substructure signif- 
icances as a function of the local cluster number density 
contrasts excluding the very extreme contrasts where sam- 
ple sizes are small and the results quite noisy. The bars 
represent the formal ler errors (excluding cosmic variance). 
It is seen that the average significances of /3 and LEE de- 



crease with density contrast. The effect is supported by 
subsamples of nearby clusters with comparatively large 
numbers of X-ray photons. In addition, we compare the 
sensitivity of the results to boundary effects possibly in- 
troduced by the galactic plane. The same trends are found 
(although with larger scatter) if only those clusters where 
the maximum distance of the 5th neighbour cluster is far 
away from any survey boundary were analyzed. 

The same effect appears also on a larger metric scale 
when the aperture radius for substructure detection is in- 
creased from IMpc to 3 Mpc (see Fig.||, right). Here we 
could follow the trend to even higher density contrasts, 
although error bars get quite large. For the latter case 
we can, however, not rule out that neighbouring clus- 
ters not necessarily in the process of merging with the 
programme cluster might artificially increase the SORs 
(notice that the average significances are in general lower 
for the 3 Mpc apertures compared to the 1 Mpc results) . 
In addition, many clusters have significance radii smaller 
than 3 Mpc especially at lower redshift so that the effect of 
background sources or local variations of the X-ray back- 
ground appears to be less clear. 

To summarize, if the (3 and LEE statistics are regarded 
as reliable indicators for substructure (see Sect.f|) then the 
results shown in Figs. @ to || clearly indicate that the frac- 
tion of substructured clusters increases with local cluster 
number density. 

Contrary to the results obtained with (3 and LEE, the 
average elongation significances are found to be almost 
insensitive to local cluster number density contrasts. A 
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Fig. 9. X-ray luminosity versus redshift for REFLEX+BCS 
clusters (small points) and for halo and relic clusters (filled 
circles) . 



decrease of (Sfel) might be 
but the effect seems to be 
value is just consistent with 
different behaviour of FEL 
However, as noted in Sect.|J 
not necessarily restricted to 
'S'fel > values are expected 
regions, which could reduce 



indicated on the 3 Mpc scale, 
not very strong. A constant 
the la confidence limit. The 
is not yet fully understood, 
the presence of elongation is 
merger events so that low < 
in both high and low-density 
a density-dependent effect. 



7. Substructures in halo, relic, and cooling flow 
clusters 

Radio halo and relic clusters are found among the most 
X-ray luminous galaxy clusters. This is shown in Fig. ^| 
where the X-ray luminosity in the energy band 0.1-2.4keV 
is plotted as a function of redshift for the REFLEX+BCS 
reference sample (dots), and the halo/relic clusters (filled 
circles). A threshold luminosity, Lx — 4.0 x 10 44 ergs _1 , 
might be introduced above which most of the halo/relic 
clusters are located. In the following we assume that the 
given halo/relic sample is representative and compare the 
observed fractions of halo/relic clusters with substruc- 
tured RASS-3 X-ray images with the corresponding results 
obtained for cooling flow and REFLEX+BCS clusters. 

For a proper comparison of the frequencies of clusters 
with substructure obtained with different samples one has 
to equalize the efficiencies of substructure detection which 
might be different for the samples. Sect.|| shows that the 
number of photons per target, Np^, is the most crucial 
factor determining the efficiency of substructure detection 
(the redshift bias will be discussed below). The first step 
thus is to equalize the frequency distributions of Af p h for 
the different samples. 
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Fig. 10. Frequency distributions of substructure significances 
of halo/relic clusters (lines) and Reflex+BCS clusters (points). 
The computation of the latter frequencies and their la stan- 
dard deviations is described in the text. 



Both the halo/relic and the cooling flow samples 
have on average higher X-ray fluxes yielding more X-ray 
photons per cluster target compared to a typical RE- 
FLEX+BCS cluster. Therefore, the REFLEX+BCS sam- 
ple is in general less sensitive to substructure detection 
compared to the other samples. The naive comparison of 
the substructure tests would thus lead to a systematic 
underestimation of the frequency of substructure for RE- 
FLEX+BCS clusters. 

In order to correct for the systematic effect one has 
to dilute a larger reference sample so that its normalized 
cumulative probability distribution function, P(< N p h), 
resembles the distribution function of a smaller test sam- 
ple. The dilution is performed as in standard Monte Carlo 
experiments by randomly selecting from a uniform distri- 
bution a number between and 1, computing the corre- 
sponding N ph value as given by the distribution function 
of the test sample, and selecting the reference cluster with 
a A^ph value next to N ph . These steps are repeated until 
the replicant of the reference sample has the same sample 
size as the test sample. In order to get more representa- 
tive results many replicants of the reference sample are 
created, giving better estimates of the mean values and 
their standard deviations. 

The resulting significance frequency distributions and 
their la errors (points with error bars, see below) are given 
in Figs. [II] to [l2|. Notice that the latter figure compares 
halo/relic clusters with cooling flow clusters. In this case 
the halo/relic sample is used as reference. The standard 
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Fig. 11. Frequency distributions of substructure significances 
of cooling flow clusters (lines) and Reflex+BCS clusters 
(points). The computation of the latter frequencies and their 
la standard deviations is described in the text. 

deviations (error bars) are computed in the following way. 
As mentioned above the reference sample is diluted so that 
its sample size and photon number distribution resembles 
either the halo/relic, or cooling flow sample. For each clus- 
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Fig. 12. Frequency distributions of substructure significances 
of cooling flow clusters (lines) and halo+relic clusters (points). 
The computation of the latter frequencies and their la stan- 
dard deviations is described in the text. 



ter type this is done 20 000 times. The standard deviation 
is estimated from the scatter of the numbers of clusters 
selected in each significance bin (a kind of Bootstrap er- 
ror). 

It should be noted that the following analysis is re- 
stricted to a discussion of the frequencies of the signifi- 
cances obtained for circular symmetry (/?), unimodality 
(LEE), and mirror symmetry (FEL) for different cluster 
types. Therefore, the introduction of a specific confidence 
level is not necessary. The results presented here do also 
not depend on the specific binning used to count the fre- 
quencies of the significance values. The following state- 
ments on the significances are based on the frequencies 
and their errors obtained in the first S bin because here 
the differences in the substructure behaviour are most in- 
teresting. 

The comparison indicates that halo/relic, cooling flow, 
and REFLEX+BCS clusters do not strongly differ in 
their substructure behaviour. The differences found in the 
RASS-3 images between these types are all on the 1-2 a 
level. However, the general picture that cooling flow clus- 
ters appear to be more regular and halo/relic clusters more 
often substructured is supported. In the following more 
specific results are described. 

Comparison of halo/relic with REFLEX+BCS clusters 
(see Fig. |kJ) : halo and relic clusters are found to have high 
frequencies for all three substructure tests below S = 0.05 
compared to REFLEX+BCS, suggesting that halo and 
relic clusters appear to be on average more substructured 
compared to REFLEX+BCS clusters. Whereas the excess 
is only marginal for f3, about la excesses are found for 
both LEE and FEL. Halo/relic clusters thus appear more 
often bi-modal and elongated compared to other cluster 
types. 

Comparison of cooling flow with REFLEX+BCS clus- 
ters (see Fig.|ll|): cooling flow clusters are found to have 
lower frequencies in the S < 0.05 bin for j3 and LEE 
whereas a small excess is found for FEL compared to the 
REFLEX+BCS reference sample. The deviations are on 
the lcr level, but in the opposite direction as found for the 
halo/relic clusters in Fig.|l^. The measurements thus sug- 
gest that clusters with large cooling flow signatures appear 
to be more often circular symmetric and unimodal. Signif- 
icant elongations of their X-ray images are found slightly 
more often compared to REFLEX+BCS clusters. 

Comparison of cooling flow with halo/relic clusters (see 
Fig. [l2]) : cooling flow clusters are found to have more of- 
ten circular symmetric and unimodal RASS-3 images com- 
pared to halo/relic clusters. Moreover, cooling flow clus- 
ters appear to be marginally more often elongated. 

8. Discussion 

A systematic study of the morphologies of RASS-3 images 
of BCS and REFLEX clusters is given. The two surveys 
provide the largest presently available X-ray cluster sam- 
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pies and are expected to yield statistically representative 
results. Our analysis shows that many observational ef- 
fects lead to systematic errors in the statistic of merger 
events. After approximate corrections are applied an es- 
timate of the SOR of 52 ± 7 percent is found. This num- 
ber might also be contaminated by chance superposition 
of point-like (background) sources. However, the effects 
are expected to be small because a comparatively small 
aperture radius of 1 Mpc for substructure measurements 
is used. 

How does this SOR estimate compares with the re- 
sults obtained with similar projects in X-rays? As men- 
tioned above, with 208 Einstein IPC images Jones & For- 
man (1999) find by visual inspection a SOR of 41 percent. 
With 65 Einstein IPC images Mohr et al. (1995) used the 
emission-weighted centroid variation for substructure de- 
tection. Kolmogorov-Smirnov tests suggests that the sam- 
ple is representative. They found a SOR of 61 percent 
for the same confidence level (99 percent) as used for the 
REFLEX+BCS sample. 

It is thus seen that the three largest presently avail- 
able systematic X-ray cluster works give SORs of about 
50 percent. However, the conservative (formal) 3<r stan- 
dard deviation of 30 percent between the three estimates 
already indicates that there is still considerable scatter 
between different samples and methods. The conservative 
interval of substructure occurrence rates 
20 < / < 80 percent (99%, confidence range) , (9) 
for nearby clusters with z < 0.15 might give a realistic 
impression of the current situation of statistical work on 
X-ray SORs. 

The next step towards a physical understanding of the 
observed SOR should be the determination of the mass 
scales of the subclumps and the dynamical time scales in- 
volved. Note that the individual contributions of major 
mergers and accretion to (^) are not given by the mea- 
surements. Obtaining quantitative estimates appears to 
be quite difficult, even if the analysis would have been 
done with better data and refined substructure tests. 

However, large sample sizes offer the possibility to cali- 
brate the substructure events at least in a statistical man- 
ner by the application of the same substructure tests to 
both observed and simulated X-ray cluster images dis- 
tributed in flux and redshift in the same way. The simula- 
tion will also include a realistic X-ray background so that 
an artifical increase of the observed SOR through back- 
ground point sources (see above) is taken into account. 
Simulations of this kind would establish the link between 
substructure as defined by the various measures and the 
dynamical state of a cluster (Schuecker et al., in prepa- 
ration). Some interesting statistical results deduced from 
the combination of observational work and numerical ex- 
periments can be found in, e.g., Mohr et al. (1995). 

Depending on the accuracy of this comparison one 
should also try to investigate redshift-dependent effects 
where no information is available yet. High-resolution N- 



body simulations of, e.g., Gottlober, Klypin & Kravtsov 
(2001) suggest an increase of major merger rates by a fac- 
tor of about 2 between redshift z = and 0.25. 

Almost independent from the observational effects 
mentioned above the relative comparison of cluster sub- 
samples show that the fraction of clusters with X-ray sub- 
structure turns out to be a function of local cluster num- 
ber density in the sense that the fraction of clusters with 
low statistical significances for circular symmetric and 
unimodal X-ray surface brightness distributions increases 
with local cluster number density. Such a substructure- 
density relation is expected when the merging of sub- 
clusters causes the observed substructures. Although the 
analyses are based on the largest presently available X- 
ray cluster catalogue, we were forced to use flux-limited 
subsamples in order to get statistically significant results. 
When larger X-ray cluster samples are available, future 
studies will allow the extraction of statistically signifi- 
cant volume-limited subsamples. This will give a better 
defined cluster number density with the additional ben- 
efit of rough mass estimates as obtained from the mass 
X-ray luminosity relation of Reiprich & Bohringer (2001). 

The substructure-density relation of clusters appears 
to be analogous to the morphology-density relation of 
galaxies. A related effect, namely that dynamically young 
optical APM clusters are more clustered than the overall 
cluster population was recently found by Plionis (2001) 
thus supporting the present findings. 

The relative abundances of radio halo, relic, and cool- 
ing flow clusters with substructures could be studied in 
more detail. Whereas halo and relic clusters tend to show 
more often substructure in RASS-3, cooling flow clusters 
show the opposite effect. Notice that the comparisons of 
cooling flow clusters either with REFLEX+BCS (Fig. |lj) 
or with halo/relic clusters (Fig. [l^) give quite similar re- 
sults supporting their robustness because almost indepen- 
dent reference samples are used. We thus regard the find- 
ing that cooling flow cluster show in general more often 
regular and unimodal X-ray surface brightness distribu- 
tions than other cluster types as quite stable. Notice that 
a fraction of 5/22 = 0.23 cooling flow clusters show either 
Slee = or/and Sp = indicating that substructured 
cooling flow clusters are to be expected. In addition, a 
fraction of 7/22 = 0.31 have <Sfel < 0-01 indicating that 
one third of the cooling flow clusters show significant elon- 
gation within the inner IMpc radius, possibly impressed 
by the cD galaxy. 

On the other hand, although halo and relic clusters 
show more often bi-modal and elongated RASS-3 surface 
brightness distributions, they appear to share more simi- 
larities with other cluster types than cooling flow clusters. 

Nevertheless, our analyses give additional support to 
the idea that cluster mergers might trigger the formation 
of radio halos and relics whereas pre-existing cooling flows 
might be disrupted by recent merger events. 
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